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1. INTRODUCTION 

The need to reduce emissions of greenhouse gases, which cause global warming and pollute the 
atmosphere, has led researchers to look for new non-polluting ways to generate energy, such as wind power. 
Within this branch of generation, one of the most widely used generators is the double-fed induction 
generator (DFIG), shown in Figure 1, which has the capacity to continue generating when disconnected from 
the main electrical grid. This alternative has grown in the last decades due to the technical advance that has 
allowed the decrease of the investment costs and the increase of the generation capacity per kWh. Although it 
must be borne in mind that the use of wind as a primary source in the generation of electrical energy can 
bring disadvantages, since the wind speed is variable every hour of the day or season of the year during the 
generation periods, this causes a limiting when connecting the turbine to the network where optimal 
performance in terms of efficiency is not achieved [1]-[4]. 

Because the power flux between the generator and the grid does not remain constant due to wind 
disturbances, which leads to variations in voltage and frequency, which are variables that must remain in fixed 
values imposed by the network to which it is connected. Therefore, in recent years we have been working with 
DFIGs, which can work at variable speeds avoiding the mechanical stress of the machine and allowing an increase 
in the energy captured; they also have a back-to-back converter that allows regulation of the voltage and frequency, 
the flow of reactive power to the network and keeps the voltage at the generator terminals constant. In comparison 
with the asynchronous generator in general, the DFIG has a more complex configuration [5]-[8]. 
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Bearing in mind that the main disadvantage of wind turbines is the fluctuation of wind speed, a control 
system is required that makes their efficiency reach an optimum point of performance. The control systems 
ensure that the operation of the system is correct as long as it works under normal conditions or disturbances 
occur, and the control implemented must guarantee greater efficiency in terms of generating electrical energy 
by discarding the speed of the turbine and adapting the speed from the rotor to variations in wind speed in 
order to achieve a greater amount of energy generated. Currently, there are different control schemes such as 
the proportional integral (PI) controller, linear quadratic gaussian (LQG) controller, robust H-infinity 
controller, discrete linear quadratic regulator (DLQR) controller, and optimal control. One of the most used 
schemes in generation systems is the PI controller where algorithms are applied that optimize the generator 
performance, although it has the disadvantage of not receiving as much energy as possible from the wind. 

In this paper, the optimal multivariable control of the DFIG is proposed, using its own structure 
assignment through the multivariable output error state space (MOESP) algorithm. The document is 
organized as: section two is about the control of multivariable systems and how to identify them from their 
inputs and outputs. Section three describes the model in the state space for the DFIG and finally section four 
shows the results obtained from the implementation of the controller by assigning its own structures. 


Wind 
turbine 


Figure 1. Wind generation system with a double-fed generator [9] 


2. METHOD 
2.1. Multivariable control 

In this subsection of the document we will talk about the characteristics of a multivariable system 
control, and how they can be identified from their inputs and outputs. Figure 2 defines the general structure 
of a controller, where the control law is defined taking into account the accumulated one for the error of the 
control signal [10]. The error of the control signal is defined as: 


e(k) = r(k) — y(k) i 
With which the extended system will be defined by: 
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Where the control signal is calculated as: 
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Here, K and K; are constant matrices, which vary according to the control algorithm applied to the system 
studied using from the (1) to the (4). 


Figure 2. Block diagram of studied structure [10] 


2.2. Assignment of structures 

In recent years the assignment of structures has been applied in different types of multivariable 
control systems, as it is a robust control technique and a high performance. Said technique can be applied 
using state feedback or feedback of the output by means of a gain. The method of allocation of structures 
basically steps: 

— Choice of a set (or sets) of possible eigenvalues (or poles) in closed-loop. 

— Calculate the subspaces of eigenvectors for the generation of eigenvectors in closed-loop. 
— Select some specific vectors of the subspaces according to design techniques. 

— Calculate the control law for the election of the own structure. 

The multivariable output error state space (MOESP) method allows determining a well-conditioned 
solution for the problem of assigning eigenvalues by state feedback. The solution obtained is such that the 
sensitivity of the eigenvalues assigned to the perturbations in the system and the gain matrices are minimized [11]. 
The objective function for optimizing the sensitivity of the eigenvalue in general can be expressed by means 
of the (5). 


a(R) 
On(R) 


Jo = nR) = (5) 


Where o is the largest singular value of the right eigenvector array of R and o, is the smallest 
singular value of the right eigenvector array of R. The formulation of the problem of the technique of 
allocation of own structures, is based on finding the law of control by feedback of states [12] (6) to (12). 


u(k) = —Kx(k) (6) 


Given the eigenvalues A = diag(A,, Ap, ..., An) and the array of eigenvectors R, there is a K such that: 


(A+ BK)R = RA (7) 
Just if: 

UT (AR — RA) = 0 (8) 
Where: 

B=[Uo UV ]lz" ol” (9) 


With U = [Uo U,] orthogonal and Z non-singular. Therefore, the profit matrix K is explicitly defined by: 

K = ZU (RAR! — A) (10) 
The law of control by feedback states can be extended by a reference as (11). 

u(k) = —Kx(k) + K;e(k) (11) 
Where K; is the steady-state gain or reference gain defined by: 

K; =C(U—A+BK)'1B (12) 


Where C(I — A + BK)~1B the matrix pseudoinverse. 


Stability by assigning structures by applying the multivariable subspace identification ... (Fernando Mesa) 


904 m) ISSN: 1693-6930 


2.3. Identification by subspaces 
The linear model of the system can be obtained by applying a subspace approach such as the orthogonal 
decomposition (ORT) method or the MOESP method of the system around an operating point [13] (13) to (18). 


x(k + 1) = Ax(k) + Bu(k) 


y(k) = Cx(k) + Du(k) (13) 


Where the input-output data is assumed to be {u(t), y(t),t = 0,1,...,N + 2k — 2} and are given with N 
large enough and k > n. Based on the input-output data [14], the input matrix passed from the Hankel blocks 
and and future entry are respectively defined as: 


u(0) u(1) -= u(N — 1) 
Te ica a T (14) 
Ge 21) HOD. 3 UO I) 
u(k) u(k +1) = u(k+N-—1) 
U; uet 1) ves 2) = u(k 4 N) (15) 
u(2k—1) a e NEED 


In the same way, the past and future of the exit is defined y,, yp. The extended observability matrix 
of order i, where i = k- 1, the inferior triangular Toeplitz matrix, and the ratio of the decomposition LQ are 
defined respectively as: 


R= (C" (CA .2c (CASE (16) 
D 0 - 0 
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Now, we show the MOESP algorithm, estimation of a multivariable system: first, to build data matrices 
Up, Uy, Yp, Yf. Second, to perform factoring LQ from the (18). Third, to perform singular value decomposition 


(SVD) to the work matrix G = [L42L43], with G = [0 J] Ñ 1 K = Û$VT. Forth, to calculate the 


estimates of A and C from (16) J; = ogt/2, Finally, fifth, to calculate the estimates of B and D from (17), 
U Lall =U H. 


2.4. Generator model 

The DFIG doubly fed induction generator is the one that is most commonly installed in many wind 
farms. In comparison with the asynchronous generator in general, the DFIG has a more complex 
configuration [15]. The DFIG is based on a principle of the induction generator and incorporates two groups 
of similar multi-phase coils of power that have independent means of excitation. The DFIG transforms the 
power of the input turbine into electrical energy. The power produced in the stator is always positive [10]. 
The Figure 3 shows the equivalent circuit of a DFIG in the reference frame that rotates at the synchronous 
angular velocity w, [16]-[19]. From Figure 3 we have the stator and rotor flows are given by (19) to (22). 


(ve = Lsls + Lml, 


Y, = Lmls + Ll, (19) 


Where L., Lr, Lm, are the self-inductances of the stator, rotor and mutual respectively. The stator and rotor 
voltages V, and V, are expressed as: 
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Vy = Rely + 2 + jo, Y, 


de (20) 

V% = R,1, + aE + j(@, — a, )¥, 

Being J, and I, the stator and rotor currents, R, and R,;are the resistances of the windings, œw, is the 

angular velocity of the rotor and w, = w1 — ws, the angular velocity of glid [19]-[21]. The DIFG uses the 

stator voltage orientation (SVO) scheme as described in the reference [22]-[25], which allows rewriting the 
voltage (21) in the rotor in the synchronous reference frame dq as: 


ee ch R 
Via = Vra + E (v, ac L Psd + or Paq) 


(21) 
1 Lm Rs 
Va = Va — a(r Peq + wr Psa) 
The system representing the state space is presented as: 
dir ioii ; 
= a Vra — Rilpg + WsoL,Lyq) 
(22) 


dIrg_ 1 1 ' 
dt oL, (Va ~ Ry lrg = ws0L,Lya) 


r[k] 


Figure 3. Equivalent circuit of the DFIG in the frame of reference with speed [19] 


3. RESULTS AND DISCUSSIONS 

For the model presented in the previous section, the outputs of the system are the currents in the rotor 
on the axis d and on the axis q, that is, a system of two inputs and one output (23) and (24). The Figure 4 shows 
the behavior of the open-loop outputs before the inputs presented in the Figure 5. 


y= |i] (23) 


Output 1 
Output 2 


Outputs (p.u.) 
Inputs (p.u.) 


Input 1 


Input 2 
4 1 L L L 4 4 L L 4 
0 1 2 3 4 5 6 z 8 9 10 
Samples 
Figure 4. System response in open loop (source: author) Figure 5. System inputs (source: author) 
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Knowing the inputs and outputs of the system, the identification of the same is applied through the 
algorithm proposed in section 2, with which the following system is obtained. 


T 2s 
m 2» 
Teo 20 
Te an 


With the previous estimated model, a controller is designed by assigning its own structures for the 
multivariable system. The following figures show the responses of the I4 and Lq streams when applying the 
driver to the system (see Figure 6 and Figure 7). Now, when analyzing the results of the output signals when 
applying the controller by assigning own structures, it can be concluded that the models respond in an 
approximate time of 0.2 s to reach the reference and a disturbance time 0.2 s. 


Reference 1 


Output 1 


Amps (p.u.) 


1 fi 1 f f 1 1 f 
0 02 04 06 08 1 {2 14 16 1.8 2 
Samples 


Figure 6. Current response J,q (source: author) 


—— Reference 2 
Output 2 


1 1 1 
0 02 04 06 08 1 1.2 1.4 16 18 2 
Samples 


Figure 7. Current response I, (source: author) 


4. CONCLUSION 

This work proposed a controller for allocation of structures for the doubly fed induction generator 
DFIG, based on the identification of the multivariable system through a MOESP subspace method, for this 
estimation it is required that the input/output data of the plant be as variable as possible and thus know the 
behavior of the plant before any internal or external disturbance, to be adopted by the identified system. 
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Once the system was identified, the controller was designed with which the results were obtained that 
allowed us to observe how the J, and Ig outputs reach the reference levels in very short times and with very 
low oscillations. Once the system was identified, the controller was designed with which the results were 
obtained that allowed us to observe how the Iq and I;,, outputs reach the reference levels in very short times 
and with very low oscillations. Therefore, it is concluded that the controller by assignment of structures is a 
good option for the control of multivariable systems. This research did not receive any specific grant from 
funding agencies in the public, commercial, or not-for-profit sectors. 
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